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Abstract 

An Ising model with local Glauber dynamics is studied under the influence 
of additional kinetic restrictions for the spin-flip rates depending on the ori- 
entation of neighboring spins. Even when the static interaction between the 
spins is completely eliminated and only an external field is taken into account 
the system offers a phase transition at a finite value of the applied field. The 
transition is realized due to a competition between the activation processes 
driven by the field and the dynamical rules for the spin-fiips. The result is 
based on a master equation approach in a quantum formulation. 
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I. INTRODUCTION 



The dynamical properties of different models as Ising model, Sherington-Kirkpatrick 
spin glasses, influence of time dependent external fields have been the subject of continuous 
interestS. Because the conventional Ising model is a static, equilibrium model, a dynamical 
generalization was first considered by Glauberi who introduced the single-spin-fiip kinetic 
Ising model for describing relaxation towards equilibrium. The dynamics is generated by a 
coupling to a heat bath at a fixed temperature. Hence the rates of the dynamical processes 
are constrained to satisfy detailed balance. Based on that the static interaction energy of 
the spins is included in the flip-rates. In this way the equilibrium phase transition inherent 
for the Ising system in (i > 1 is also manifested in the dynamical approach such as within 
the behavior of the relaxation time and the correlation function. 

In the present paper we discuss the possibility that the kinetic Ising system undergoes a 
phase transition although the mutual static interaction between the spins will be neglected 
completely. The transition is purely mediated by the underlying kinetics. Special attention 
will be paid to dynamical restrictions for the spin-flip rates instead of including an inter- 
action. As a consequence of such constraints combined with the orientation of the applied 
field (activation energy) there occurs a competing situation leading to the aforementioned 
kinetic induced phase transition. 

To be specific let us consider a spin system in an external time independent field which 
is subject to a dynamical process (spin-flip process). However, it is assumed that the flip 
rates does not depend on the neighboring spin configuration via the static energy between 
adjacent spins manifested by a Hamiltonian. Alternatively, kinetic restrictions are included 
independently on their energy. The environment of a fixed spin is oriented in a self-adaptive 
manner only by the kind of kinetic processes. Formally, this can be achieved by chosing 
the flip rates depending on the state of the neighboring spins. Furthermore, the flip rates 
are also determinated by the temperature and the activation field. To formulate the master 
equation under the influence of constraints we use a representation of the underlying kinetic 
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equation in terms of quantum Pauli-operators. Based on the analogy to quantum mechan- 
ics the problem can be studied in a compact forml, for a recent review seei'l. The method 
allows the inclusion of the above mentioned restrictive conditions in a transparent manner 
which had been demonstrated for the so called Fredrickson Andersen mo delii, seeu. 



II. THE MODEL 

Let us use a lattice gas description with the occupation number at a certain lattice site 
defined by rii = 0, 1. Due to the relation to the spin variable, Si = 1 — 2nj, an empty 
site corresponds to an up orientation whereas an occupied site is related to spin-down. 
The local flip dynamics can be characterized by an evolution operator written in terms of 
Pauli-operator^ 

Li = A(l - dl)d, + 7(1 - di)4 (1) 

The first term gives a nonzero contribution in case of a flip from a spin-down to a spin- 
up state (or from an occupied to an empty site) with the rate A whereas the second term 
including the rate 7 is responsible for the inverse process, i.e. the flip from spin-up to spin- 
down states. Obviously, the rates A and 7 can depend on the temperature of a heat bath 
coupled to the system, on the mutual interaction between the spin degrees of freedom, as well 
as on the neighboring spin configuration in a kinetic manner discussed below. The coupling 
to a heat bath is easily realized by the following modification of the evolution operator! 

L = uJ2[{^- di) exp{-(3H/2)dl exp{(3H/2) + (1 - dj) exp{-(3H/2)di exp{(3H/2)] (2) 

The procedure is similar to a 'thermalized' Heisenberg picture in quantum mechanics. As 
the simplest case let us consider a coupling to an activation field of the strength e 

H = Y.eni (3) 

i 

If £ > for each lattice site the spin-up orientation is preferred. In the opposite case 
£ < the spin-down state is more favored. The algebraic properties of Pauli-operators can 



be used to calculate for instance exp {-pH/2)dlexp{(3H/2) = 4 exp{l/2(3e). Comparing 
eq.(|l]) and eq.@ the thermalization leads to temperature dependent activation rates which 
can be written in the case considered by eq.(0) as 

X = uexp{6/2T) 7 = z/exp(-e/2r). (4) 

It should be remarked that the transition rate from down to up state is realized at low 
temperatures with a larger probability than the opposite process (in case e > 0). For high 
temperatures, T — oo, both flip-processes appear with the same probability. Before the 
model is extended by introducing kinetic restrictions let us demonstrate its applicability. 
The evolution operator L obeys the equation 

d, I F{t)) = L I F{t)) (5) 

where the state vector | F{t)) is related to the probability distribution function P{n,t) 
according to P(n, t) = {n\ F{t)) with a complete set of basic- vectors | n) composed of Pauli- 
operators. The probability distribution function is assumed to follow a master equation 
written in the symbolic form 

dtP{n,t) = LP{n,t). (6) 

The evolution operator L in eq.(^ can be mapped onto the operator L in the mas- 
ter equation. Such an approach had been introduced by Doii using Bose operators^. 
A generalization to Pauli-operators had been proposedS^, for a recent review seeSi. 
FoUowingiSllUll'El the probability distribution P{n,t) is related to a state vector | F{t)) 
in a Fock-space according to P{n,t) = {n \ F{t)) with the basis-vectors | n) composed of 
second quantized operators. A further extension to an p-fold occupation number is possible 
and had been discussed in0. 

The relation between the quantum-like formalism and the probabilistic approach is given 
by 

\F{t)) = j:P{n,t)\n). (7) 



As it was shown firstly by Doiu the average of an arbitrary physical quantity B{n) is defined 
by the average of the corresponding operator B{t) 

{Bit)) = J2Pin,t)Bin) = {s \ B \ Fit)) (8) 

rii 

with the state function (s |= I- Using the relation {s \ L = the evolution equation for 
an operator A can be written 

dt{B) = {s\[B,L]\Fit)). (9) 

It should be stressed that all the dynamical equations covering the classical problem are 
determined by the commutation rules of the underlying operators and the structure of the 
evolution operator L. In our case the dynamics of the model is given by spin-fiip processes 
indicating a change of the local occupation number or the spin orientation, respectively. 
Based on the above introduced approach the evolution equation for the occupation number 
operator or the averaged spin for the simple fiip process eq.(|^) is easily derived. We find 

u'^d{n) = -(n)e^/2T ^ _ n)e-'/^^. (10) 

The solution in a spin representation is well known and reads 

(^(t))=tanh(^) + S'oexp(-t/r) with r = ^^J^. (11) 

The first term represents the stationary state whereas the second one describes the exponen- 
tial decay with a relaxation time r. For low temperatures one finds r oc exp(— ^). In the 
opposite limit the relaxation time shows an algebraic dependence on the activation energy 
T OC |(1 — 2(|;)^). The results are independent on the sign of the activation energy. 
Now we discuss an extension of the model introducing constraints within dynamical rules 
specified below. It will be demonstrated that in this case the sign of e leads to a complete 
different behavior of the system. To this aim an anisotropic model is studied where the 
kind of anisotropy is not originated by the static interaction. Instead of that the fiip pro- 
cesses will be modified in the following manner. Whereas the fiip process from spin-down 
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to spin-up orientation should be possible whenever a double occupancy is avoided which 
is automatically guaranteed by the algebraic properties of the Pauli-operators, the inverse 
process, spin-up to spin-down flip, should be restricted. Such a flip is only allowed when the 
neighboring configuration fulfills the condition 

^E(l + ^.)</ (12) 

A local spin flip at the lattice site i is only permitted when the number of nearest neighbors 
denoted by satisfies the last condition. Here / is a certain but fixed number specified 
below. Physically the confinement means that flips from up- to down-spin orientations are 
reduced when the total number of adjacent spins with an spin-up alignment does not exceed 
the number /. Such kinetic restrictions are already introduced by Fredrickson and Andersenil 
within a f-spin facilitated kinetic Ising model which is a candidate describing the ultraslow 
dynamics in the vicinity of a the glass transitionlii"0i. 

Here, we analyse an anisotropic version of the model by assuming that only one of the flip 
rates, A or 7 are modified which is different from that one discussed for glassy materials. 
The kinetic restriction manifested by eq.(|r^) means that a local change from a spin- up state 
to a spin-down state is not possible when more than f nearest neighbors with respect to a 
given site are in the up position. Mathematically, the confinement is realized by assuming 
that the corresponding transition rate, here 7, will be replaced by 7 — > 771^^ . . . rir^. As the 
simplest choice let us chose f = z/2 where z is the coordination number. This leads to the 
following evolution operator 

L = ^^Xih-jf (1 - di) exp{-PH/2)dlexp{PH/2)nj^ ■ ■■rij^ 

+ [(1 - 4) exp{-f3H/2)di exp{f3H/2)] (13) 

where Xiji-jf is only nonzero when all the indices indicates nearest neighbor sites. Here 
we consider the case f = 2 which corresponds to a cubic lattice in two dimensions. The 
general results seem to be independent on the special realization. As the main feature of 
the model there occurs a conflicting situation between the static preferred orientation of the 
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spins manifested in the sign of the activation energy e and the kinetic rules. That case will 
be discussed in detail in the forthcoming section. 

III. ANISOTROPIC MODEL 

Using the algebraic properties of Pauli-operators the evolution equation for the averaged 
occupation number reads 

dt{nj) = uJ2Xjrs[ewi-£r/'2T){{l - nj)nrns) - exp(£j/2r)(n^)] (14) 

r,s 

As already demonstrated i nl the mesoscopic dynamic equations within the classification 
introduced by Hohenberg and Halperinill is rederived in a coarse grained approximation. In 
the present case the equation is written in the form 

dt{n^) = -e^l^^z{z - l)(n,-) + e'^/^^ ^ K) - (15) 

rs 

It consists of a local part (linear term) and nonlocal terms. After Fourier transformation 
and introducing the abbreviation {n{k,t)) = nk{t) the last equation reads 

v~^dnk{t) = -e''/'^nk{t)z{z - 1) 

+ e"^'''^ [-tpqnk-p^q{t)np{t)n^{t) + efc_ppnfe_p(t)np(t)] (16) 

Here we have introduced the dispersion relation (/ is the lattice constant) 

Cpq = 4[cospa;/ COS Qyl + cospyl COS q^l] + 2[cos(pa; — qx)l + cos(py — qy)l] (17) 

The homogeneous stationary state uq obeys the equation 

= -Xno + 7(^0 - ^o) (18) 
This equation yields three stationary points 

^ = n2/3 = ^[l± Vl-4p] with p = ^ = exp(5/T). (19) 
Making the conventional ansatz n{k,t) = nQ6{k) + m{k,t) we get 



— * — * — # 

dtm{k,t) = —A{k)m{k,t) with 

A{k, t) = Aeoo + 7[(2efco + eoo)^o " 2efco%] 



(20) 



From the hnear stabihty analysis we conclude that the spin-up orientation rii = is always 
stable independently on the sign of the activation energy e. For the nonzero solutions the 



Obviously, the sign of gap A(0) in eq.(^) is responsible for the stability. Whenever e > 
the both nonzero solutions n2 and are unstable. Hence, we find that in this kind of an 
anisotropic model the up orientated spin state is the single stable one. The result is expected 
because the static field (activation energy) favors the spin-up state, e > 0, and the dynamic 
rules support the same orientation with a larger probability as the opposite one. The result 
is independently on the temperature and the strength of the activation energy. 
The situation is complete different analysing the case e < 0. The static energy would 
favor the down orientation of the spins due to the sign of the field. In contrast to that the 
dynamical rules tend to avoid such a state because states with a down-spin orientation are 
reduced by the kinetic restrictions. As the result of the competition between both processes 
there appears a critical temperature Tc= -j^. Above this temperature only the state ni = 
is a stable one. For T <Tc both states rii and n2 are stable, where n2 can be rewritten 



The situation described in the case e < is remarkable. For high temperature the spin- 
up state is realized. Whereas the flip rates in this case are of the same order the static 
term would favor a spin-down state. Due to the kinetic restrictions such a state is not 
realized, or with other words the dynamic rules with the above introduced constraints work 
much stronger than the activation field. The phase transition is exclusively induced by the 



function A{k) introduced in eq.(pO|) can be written in the form 



A{k) = A(0) + Aeoo(2 - cos{kJ) + cos{kyl)) 



with A(0) = €007(^^0 - 2p) 



(21) 




(22) 
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kinetics. For very low temperatures the flip-rate to an up state (promoted by A) is very small 
compared with the rate for the opposite process. Insofar, the kinetics favor the occurrence 
of the spin-down orientation. From this point of view one should expect that for T < Tc 
only the state ^2 7^ is realized. This intuitive result can be confirmed by comparing the 
energy of both states. Starting from eq.(|T^) we find an energy function F which reads 



—n n H — n (23) 

where we have considered only the homogeneous case which seems to be sufficiently. The 
function F plays the role of the free energy. However, it should be stressed that F is 
not simply a Ginzburg-Landau functional because different to that one the prefactor of 
the quadratic part does not change its sign at the phase transition temperature (A > 0) 
reflecting the fact that the system does not possess a phase transition in equilibrium. Now 
let us calculate the energy difference AF = F{n2) — F{n = 0) with the result 

AF < for T < To 



ln(4.5) 

AF>0 for T>To (24) 

From here we conclude, that the solution ni = 0, spin-up orientation, is energetically pre- 
ferred within the interval Tq < T < Tc whereas for T < Tq the nonzero solution n2 is stable 
and is related to the lowest energy. The situation is comparable with a flrst order phase 
transition, but as already stressed, the transition is originated by the kinetic restrictions and 
it is not driven by the energy. 



The above introduced model can be easily extended by considering other restrictions. For 
instance the flip process spin-down to spin-up can be reduced when / = 2 neighbors are 
in the down position. In our approximation it means that the transition rate A is replaced 
by A(l — ny. The result does not change when simultaneously one substitutes A by 7 
and 6 by —e. Another extension can be realized when both flip processes are conflned 
for instance by assuming A A(l — n)^ and 7 7n^. When the activation energy is 
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positive, £ > 0, one finds three stable equilibrium solutions rii — 0,n2 — 1,^3 = with 
F{ni) < F{n2) < F^n^i). In the opposite case of a negative activation field £ < we get two 
stable solutions ni = and n2 = 1 with F{n2) < F{ni). A third solution = is only 
stable for T > Tc- Moreover, other reahzation for constraints does not lead to significant 
changes. 



IV. CONCLUSIONS 

In the present paper we have considered a system with competing forms of interactions. 
Whereas the static energy is simply the energy for spins in a field which does not lead to 
an equilibrium phase transition the dynamical rules for the thermal activated flip processes 
are supplemented by additional restrictions. These restrictions are self organized by the 
orientation of neighboring spins with respect to a given one. Whenever the field and the 
kinematics favor the same spin orientation a phase transition is suppressed. In the opposite 
case, when orientation preferred by the static activation field will be restricted or even when 
a certain order is prevented by the dynamics the system offers a phase transition. There 
appear two stable solutions below a critical temperature Tc where this temperature is purely 
determined by the field and not by quantities included in an energy functional of Ginzburg- 
Landau type. To be specific, we have analysed the case that a down-orientation (occupation 
number n = 0) is supported by the orientation of the field e. The fiip-rate A responsibly for a 
spin-fiip from down-to up-orientation is purely given by the activation process. However the 
inverse process, a spin-fiip from up- to down-orientation (rate 7) is additional constrained 
by the dynamical alignment of neighboring spins. In particular, the mentioned fiip process 
is strongly restricted when more than a half of the neighboring spins of a fixed one is already 
in an up-position. This confiicting situation leads to the phase transition at Tc. For low 
temperatures the relation 7 A is valid. The flip process, spin-up to spin-down, is realized 
with a much higher rate as the inverse process. Because the restrictions favors an up- 
orientation we find for a finite temperature that the state n ~ 1 is the stable one. This state 
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is always realized for T < Tq. However, there exists a second stable solution with spin-up 
orientation but with a higher energy. On that manner this state is only a metastable one. 
At To both states are simultaneously stable and coexist with the same energy. For T > Tq 
the up-orientation is favored with a lower energy. When T > Tc only the up-orientation is 
stable, that means the restrictive dynamics determines the behavior of the model. 
It seems to be for the first time that such a kinetic mediated phase transition is discussed. 
As a further step we should take into account fluctuations to observe the detailed behavior 
in the vicinity of the branching point. 
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